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Explicit Gauge Fixing for Degenerate Multiplets: 
A Generic Setup for Topological Orders 
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We supply basic tools for the study of the topological order of a multiplet which is an 
eigenspace of a finite-dimensional normal operator with continuous parameters. We allow in- 
trinsic degeneracies within the multiplet where a well-known standard procedure does not work. 
As an important example, we give novel expressions for a spin Hall conductance for unitary su- 
perconductors with equal spin pairing. Generic topological orders will be treated in this unified 
manner particularly with nontrivial topological degeneracies. 
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It has been gradually clarified that many physically 
important phenomena have origins in their topological 
orders. 1 ' 2 Some of them include (fractional and inte- 
ger) quantum Hall effects, 3-12 Haldane spin chains, 13 
solitons in polyacetylens, 14 anisotropic superconductors 
and superfluids, 15-20 chirality order in an itinerant mag- 
netism, 21 spin transport (spintronics) as a realistic ap- 
plication of Thouless pumping, 22,23 polarizations in in- 
sulators, 24 and the exotic electronic states of graphite. 25 
Strong correlations between electrons cause exotic mean 
field states 26,27 and effective quasiparticles such as com- 
posite fermions 28 which can also be discussed in terms of 
the topological orders. In many cases, the topological or- 
der itself is hidden in bulk systems but exhibits apparent 
physical consequences at the boundaries of the systems, 
such as in edge states of the quantum Hall effects, 6,29,30 
local moments near impurities in the Haldane spin chains 
(Kennedy triplets), 31 vortices and zero- bias conductance 
peaks in anisotropic superconductivities and boundary 
local moments in carbon nanotubes. 25 

In many cases, nontrivial topological orders appear by 
restricting their physical space in a manner in which a 
type of gauge structure naturally emerges. 32-34 To char- 
acterize the quantum state of a specific system, one must 
explicitly determine gauge invariant quantities for the 
physical states. The (first) Chern number 8,35,36 is such 
a candidate and it has been used for several characteri- 
zations of topologically nontrivial states. 7,30,37,38 

In this paper, we present a generic setup for the dis- 
cussion of the topological order explicitly, particularly 
focusing on gauge fixing. A standard procedure for fix- 
ing gauge was reported by Kohmoto. 8 This is well known 
today. However, the procedure does not work when de- 
generacy exists. If degeneracy is accidental, that is, it 
exists at certain special parameter values, it is negligi- 
ble. However, in several interesting situations such as in 
unitary superconductors, degeneracy is due to an intrin- 
sic symmetry, that is, the standard procedure cannot be 
applied for any values of the parameters (see below). In 
such cases, the present generic gauge fixing procedure is 
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essentially important, particularly for numerical calcula- 
tions of Chern numbers. We extend the standard proce- 
dure to general situations which allow intrinsic degenera- 
cies of eigenstates. A typical situation where our method 
is crucial is the calculation of spin Hall conductances 
for numerically obtained BCS Hamiltonians, where the 
order parameters are given numerically by minimizing 
the mean field free energy. Then quasiparticle states are 
obtained by diagonalizing a Bogoliuvov-de Gennes equa- 
tion. When the order is unitary, it has an intrinsic de- 
generacy 39 which prevents direct applications of the stan- 
dard procedure to the calculation of the spin Hall conduc- 
tance. Also, when a physical ground state has a nontrivial 
topological order and it lives on a genus g{> 0) Riemann 
surface, a fundamental topological degeneracy can occur 
with degeneracies q 9 with some integer q. 2 A typical situ- 
ation is the fractional quantum Hall effect with the filling 
factor v = l/<7- 40 In such a degenerate case, the present 
extension is indispensable. Further generic expressions in 
the present paper can be applicable to a wide range of 
physically interesting situations. 

Multiplet and Unitary Equivalence: Let us consider 
taking a normal operator L, {Lr L — Lit'), in an N- 
dimensional (N < oo) linear space. This implies that L is 
diagonalizable by a unitary matrix, U = {ipi, • • • , ^n), as 
LU = U£, £ = diag (ei, • • • , ejv)- Note that normal op- 
erators include hermite, skew-hermite, unitary and skew- 
unitary operators. (Correspondingly, the eigenvalues ej's 
are real, pure imaginary, and on the unit circle on the 
complex plane). Also, we assume that the operator L 
is labeled by a set of continuous parameters as L(x), 
x = (xi, • • • , Xd=dim v)j where V is a ri-dimensional pa- 
rameter space. 41 Various physical realizations of the op- 
erator L are (i) momentum-dependent Hamiltonians in 
the quantum Hall effect 7,8 and an anisotropic super- 
conductivity, 17,42 (ii) parameter-dependent Hamiltoni- 
ans L(x) = H(x) in the discussion of Thouless pump- 
ing 22 and the Berry phase, 33,34 and (iii) a time evolution 
operator, L(x) = Texp(-(i/fi) / dtH{t)). 

Now construct an M-dimensional multiplet (a linear 
space) W(x) with the parameters which we considered 
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(dimVK = M < N). Take M linearly independent or- 
thonormalized bases tpi, i — 1, • • • , M for W{x) (M col- 
umn vectors of dimension N) and form aniVxM matrix 
(the basis of the multiplet W) as 

*(x) =0>i,ip2,-- ■ ,i>M)- 

Supplementing the basis of the orthogonal multiplet 
Wq (W FFc = K w ), they form an orthonormalized 
complete basis for R w as (*, * c ) = 

Here, we must be cautious. We require the multiplet 
W{x) to be uniquely specified by the parameter x. How- 
ever, this does not necessarily mean that the basis H>(ar) 
is uniquely specified by x, which gives us the freedom to 
change the basis. This leads to ambiguity in specifying 
the basis of W(x). That is, we may take a different ba- 
sis, \I>'(x) = \l>(x)u>(x), where u> is an M-dimcnsional 
unitary matrix. The ambiguity of the gauge freedom u> 
can be clarified, for example, when one tries to construct 
the multiplet from the ipi's obtained numerically. This 
provides a gauge freedom for a connection we shall de- 
fine. This was first observed by Wilczek-Zee in a study 
on non-Abelian Berry phases. 34 

To specify the multiplet W(x) uniquely using x, the 
multiplet should include all of the degenerate eigen 
spaces as 

€ t (x) ^ €j(x), 

i € I in = {1, • • • , M }, j € I out = {M + 1, ■ ■ ■ , N}. (*) 

We describe this condition the existence of a generic en- 
ergy gap. The degeneracy within the multiplet is allowed, 
ti = Cj,i,j G /*", which exists, for example, in unitary 
superconductivity 42 and spin degenerate cases. The de- 
generacy between the multiplet W we considered and the 
supplementary Wc is not allowed. 

The parameter space and the multiplet are specified by 
the concrete topological orders we shall study. Let us list 
some examples of Vs with corresponding multiplets Ws: 
(1) V: the Brillouin zone, (1-i) W: a collection of Lan- 
dau level wave functions in the quantum Hall effects, 7 
(1-ii) W: quasi-particle states in anisotropic supercon- 
ductors with and without equal spin pairing, 17 and (2) 
V: a set of external parameters in the study of Thouless 
pumping 22 and the Berry phase, 33 ' 34 W: generically de- 
generate ground states, to be specific, for example, V: 
(2-i) a collection of fluxes passing through the systems, 
Aharonov-Bohm fluxes and the strength of vortices in the 
type II superconductivity and (2-ii) parameters specify- 
ing the axes of spins in magnetic systems. 

Connection and First Chcrn Number : Define a non- 
Abelian connection one-form A which is an M x M ma- 
trix as A = ^d^. Correspondingly, A' = = 
u)~ 1 Au) + uj~ 1 duj. This is a gauge transformation in our 
problem. Also, define a field strength two-form = 
dA + A 2 which is transformed by gauge transforma- 
tion as J 7 ' = u^Tu. Then Tr^" is gauge invariant 
and an integral of Tr J- over the two-dimensional ori- 
entable compact surface S (dS — 0) in V gives the first 
Chern number 35 ' 36 C s = ^ f s TrF = ^-J s TvdA. 
For example, consider the three-dimensional Brillouin 
zone V = T 3 (k X7 k V7 k z ) as the full parameter space and 



the two-dimensional Brillouin zone, S = T^ y (k z ), with 
the fixed third momentum k z for the surface S. 43-45 

A global connection over the full surface S is not 
allowed to exist in a system with a nontrivial topo- 
logical order. Then let us divide the integral region S 
into several patches Sr (S = L)Sr, R = 0, 1,2, •••), 
where the connection An is locally defined within 
Sr as A(x) = ^4^(a;). 8 ' 35 Furthermore, we assume 
each Sr,R — 1,2, does not share any bound- 
aries with dSo = — Sij>i 9Sr. When the connection 
Ar is related to A® by the gauge transformation u>or 
as Ar = oJq R Aq(jJor + u}q R (1u:or, the Chern num- 
ber Cs is written using the Stokes theorem as Cs = 
57 Eii>i Io Sr ImTrw^dwoij. 

Explicit Gauge Fixing: Topological invariants are usu- 
ally given by gauge-dependent quantities. To evaluate 
the expression, one must fix the gauge. Without fixing 
it, wc cannot have any well-defined derivative. Now let 
us explicitly fix the gauge for the multiplet. Although 
the basis \I> has a gauge freedom, a projection opera- 
tor into the multiplet P = is a gauge invariant. 
Define an unnormalized basis, from a generic ba- 
sis, * (an N x M matrix ), as = P* = *77 $ 
(t7 $ = \J r ' i '<l>). The overlap matrix of the basis 0$ = 

^£^1% is generically semipositive definite. Then, only 
if the determinant of the matrix 0$ is nonzero, we 
can define a normalized wavefunction, = Vl/^o^ 1 , 

where o$ = u\ , diag(\/Ai, ■ ■ ■ , V%m)U$ with 0$ = 
C/$diag(Ai, • • • , \m)U&. This o$ is hcrmite and posi- 
tive definite. Now we define the connection A§ with the 
gauge fixing by 4? as A<s, = ty^dty®. This is well defined 
unless dete>$ = 0. 8 - 37 

Define regions S R ,R = 1,2, as (infinitesimally) 
small neighborhoods of zeros xfj's of detC$(x) and S 1 * 
as a rest of S as 

S=\JS% t detO^){^ at^ e 5l?/il,2,.... 

We use this gauge for the region Sq, and for the region 
Sft, we use a different gauge by <I>, with detO|, ^ 
everywhere in S R . The transformation matrix between 
and is obtained as u> — o^rj^rj^o^ 1 . Since o$ 
and 0|, are strictly positive definite at the boundaries 

8Sq, we have ImTr logo; = — ImTr \og& P<&. Finally, 
we obtain an expression for the first Chern number with 
explicit gauge fixing as 

C s =-NT(S) = -J2nE(S%) 

R>1 

nn(S%)=—(f d'Sl, ft = *) = Arg det * f P*, 

where Nq(S) is the total number of signed vortices with 
the vorticity nn(Sft) inside the region Sr, R > 1 ( 
<9S* = - Ur>x dS%). i6 Since Q = Arg det »/t det rj^ , 
all the vortices of fi(<&, <&) are given by zeros of det^, 
(x*) and detr7 $ , (x*). The Chern number is obtained 
by summing up the vorticity only at x*. This form of 
Chern number is not found in literature. 
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Since we assume that the two-dimensional surface S is 
compact and Q is regular except at xf , • • • and xf , ■ ■ ■ , 
a union of curves (UrSr) U (L)rS%) is contractible to 
a point within a region where is well defined. This 
implies that N T - T (S) = N T T - (S). That is, the 

vector field O depends on the gauge (choice of «& and <&) 
but the total vorticity Nq(S) is a gauge invariant of the 
multiplet W. None of the vortices has any direct physical 
meaning. Only the total number of vortices Nq(S) has 
a physical significance. 

The projection P is essential for carrying out the 
present gauge fixing procedure. It has also an integral 
representation, P = §yGt{z), where Gt = (In — 
L)^ 1 and the closed curve T encloses all of the eigen- 
values ei, • • • , ejif inside, but not those em+i, ■ ■ ■ ,£jv on 
the complex plane. 47 From this form of projection, the 
stability condition of the generic gap (*) for obtaining 
a well-defined multiplet is clear. The first Chern num- 
ber has an apparently gauge-independent form given by 
Tr:F = -Tr dPPdP as well. 48 Also, the Chern number 
for the multiplet is expressed as 



Cc 



=-£ [< 



dtf{G c (e k )}*{G c (e k )}dL) 



4>k i 1 



where Gc = (I — P)Gt- This is equivalent to a Kubo 
formula in the case of the quantum Hall effect. 7,48 This 
formula is particularly important since mathematical ob- 
jects such as Chern numbers have a direct relation with 
a physical quantity such as a Hall conductance. Surpris- 
ingly, this is observable in a bulk system. 

Sum Rules: Assume the multiplet W is a direct sum 
of orthogonal multiplets W\ and W 2 as W — W\ © W 2 , 
which is expressed by bases and ^2 ( orthonor- 
malized in each multiplet) as *!? — (\I/i,\I r 2 ), where 
= I Ml , *2*2 = Im 2 , *I*2 = Ml M 2 , and 

>M 2 Mi ■ 



The connection is given as A ■ 



(d*i,d* 2 )- 



Thus, a trace of the connection is additive as Tr A 



TrAi + Tr A 2 , where Ai = and A2 = *2 rf *2- 

From this simple observation in the connection level, a 
sum rule for Chern numbers is as follows: Cs (W\®W 2 ) — 
C'siWx) + C S (W 2 ). The sum rule in the field strength 
level was previously discussed. 49 A simple consequence 
of the present sum rule is the total sum rule, that is, 
the Chern number of the total multiplet Wt always van- 
ishes, 42 E^Wi) = C s (®iWi) = 0, since P^ Wi = 
In- 

One-Dimensional Example (dim W = 1): When the 
multiplet is one-dimensional, such as ^ — ip and (P)ij = 

ipiipj, we have <& P& = (P)ni = ip* N ipi by taking * 
and * as = (1, 0, • • -0) and = (0, • • •, 0, 1). Then 
the Chern number is given as 



C s = - — f dim log^i /ip N ) 
27r Jds* 



ers two-dimensional electrons on a lattice with the flux 
<j> per plaquette, one-particle states are given by q bands 
when <j> = p/q with the mutually prime p and q. Further- 
more, the spectrum is given by the famous Hofstadter's 
butterfly. When the fermi energy Ep is in the j-th en- 
ergy gap from below, the Hall conductance a is given by 
the sum of the Chern numbers of the j bands. 7 In this 
case, take a multiplet from a filled fermi sea (W = FS) 
and construct the basis of the multiplet from the j Bloch 
states i>j(k) below E F as \& = (if>i, ■ ■ ■ , ipj), M = j, then 
the Chern number C^2 (FS) naturally gives the Hall con- 
ductance (j xy which is the sum of the Chern numbers of 
the filled bands. 7 ' 37 - 38 ' 49 

Dirac Monopolc: When the dimension of the total 
Hilbcrt space N is 2, only the nontrivial multiplet is 
one-dimensional M = 1. Then take an hermite Hamil- 
tonian H(x) = R(x) ■ er for the normal operator L 
where er's are Pauli matrices and R(x) is a real three- 
dimensional vector ( (R, 9, <f>) is a polar coordinate of 
R). As an example, consider the multiplet with the 
energy -R as 33 **_(a;) - (- sin e**™ cos ^). 
The projection is given as P = * -M/l = 



where S% includes a single zero of detO$ = l^il 2 - This 
is a well-known classic expression. 8 

Multiplet of Several Landau Levels: When one consid- 



-e"^ sin f cos § 



cos 



.2 8 



. Using a gauge 



(cos \ , sin J ) and * <& 



secular equation, Hip = Etjj, H = 



sin 2 § 
-e^ sin f cos § 
by <I> and l> as *<!> 

(cos |, sin |), dct C$ = and detOj, = give 

(6(x),4>(x)) = (XiO arL( i (X)0i respectively. This 
clearly shows that the positions of the vortices de- 
fined by the vector field = Arg(— sin | cos ^ + 

e+'McosfsinfX-sinf cosf + e -*W-« cos f sin |) 
are gauge-dependent and do not have any direct phys- 
ical meaning. One can chose the positions of the vortices 
as one wishes. 

Unitary Superconductors: Let us first consider the 
simplest case, that is, the unit cell includes only one site. 
Then the Bogoliuvov-de Gennes equation for generic su- 
perconductivity is given in a momentum space by a 4 x 4 

el 2 A 
At -eJ 2 

As for the unitary order, the order parameter matrix A 
is written as A = |A|A , |A| > 0, AA f = |A| 2 / 2 , 
where Ao is a 2x2 unitary matrix. Then the eigen- 
states (quasiparticle) are doubly degenerate as t ip-(w) — 
(— sin I w, cos I A «;, ), for example, for the E = —R 
state where w is a normalized arbitrary two-component 
vector, w^w = 1, R = \/\A\ 2 + e 2 , e = Rcos9 
and |A| = RsinO. Now let us construct a multiplet 
for the degenerate E = —R quasiparticle bands as 
1 3/_ = ((/)_(idi), i/)_(m)j)), where w\ and W2 form an ar- 
bitrary two-dimensional orthonormalizcd complete set: 
w\wj — Sij, W\w\ + w 2 «4 = I 2 - Then the projec- 
tions P - are given in a gauge invariant form as P = 

J 2 sin 2 | -Aosinfcosf \ 
A + . fl fl r o a ■ Now let us fix 

A^ sin \ cos I I2 cos 2 | J 

the gauge by choosing =' (0 2 , 1 2) and *€> =' (J 2 , 2 ). 
Then we have detC_ t $ = cos 4 |, O = ArgdetA = 
Argdet A. Since the overlap determinant detO-,$ = 
vanishes at 9 = it for the multiplet W- , the Chern num- 
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ber is given by 

Co = - — ^ / d'Argdet A, 

where p's are points on the surface S which are specified 
by 6 = 7T , that is, |A|(p) = and e{p) = -R{p) = -E(p). 
This is a novel expression for the generic spin Hall con- 
ductance for the unitary superconductors with equal spin 
pairing. In the previous work, the Chern number was 
given as the sum of two integers using an eigenvalue 
equation for the unitary matrix Ao. Here, we give a 
direct expression only using the order parameter ma- 
trix A. Furthermore, if one parameterizes the unitary 
2x2 matrix A as A = e je e m ' <T , |n| = 1, we have 
Cg = — i J2 P § dQ. The present method is crucially im- 
portant and efficient when the order parameter is given 
by numerically solving a BCS self-consistent equation 
with a large unit cell. Even in this generic situation, to 
evaluate the spin Hall conductance, we must determine 
arbitrary orthonormalized (degenerate) eigenstates by 
diagonalizing the mean field Hamiltonian (Bogoliuvov- 
de Gennes equation). 

Finally, we mention the higher order Chern numbers 
C n {S n ) = N n J sn Tr^ n ,n = 2,3,- ••, where N n is a 
normalization constant and S n = S x • • • S, (n- times). In 
principle, they can also explicitly be evaluated using the 
present gauge fixing procedure as a sum of integrals over 
the (2n — l)-dimensional spheres S 1 * enclosing (2n — 2)- 
dimensional regions P* which are defined by the zero of 
det 0$ in S n . They should also help in the characteriza- 
tion of the topological order in complex situations. 
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